THE FUNDAMENTAL GAP CONJECTURE FOR POLYGONAL 

DOMAINS 



ZHIQIN LU AND JULIE ROWLETT 



Abstract. In 1985, S. T. Yau made the following "fundamental gap conjec- 
ture," |25l . For a convex domain Q C M", 

(0.1) 5(n) := (A2(0) - Ai(n)) > 371^ 

where d is the diameter of the domain, and < Ai(f2) < A2(f2) are the 
first two eigenvalues of the Euclidean Laplacian on f2 with Dirichlet boundary 
condition. The scalar invariant ^ is the gap function. We restrict attention to 
planar domains. Our main result is a compactness theorem for the gap function 
when the domain is a triangle in M? . This result shows that for any triangles 
which collapse to the unit interval, the gap function is unbounded. Due to 
numerical methods, we expect that the fundamental gap conjecture holds for 
all triangular domains in R-^. We show with examples that the behavior of 
the gap for collapsing polygonal domains is quite delicate. These examples 
motivate a technical result for collapsing polygonal domains giving conditions 
under which the gap function either remains bounded or becomes infinite. 
Our work initiates a general program to prove the fundamental gap conjecture 
using convex polygonal domains. 



1. Motivation and results 

For a Schrodinger operator on a compact convex domain, after the first eigen- 
value the next natural object to study is the gap between the first two eigenvalues, 
known as the fundamental gap. This includes the work of [25], [T3], [TD], [5], [2S], 
EZ]) III]> and many other authors. While it is always interesting to understand 
the interaction between the eigenvalues of a differential operator and the geometry 
of the domain, beyond purely mathematical implications the fundamental gap has 
physical implications. For the heat equation, the gap controls the rate of collapse 
of any initial state toward a state dominated by the first eigenvalue and is of cen- 
tral interest in statistical mechanics and quantum field theory. In analysis, the 
gap is important to refinements of the Poincare inequality and a priori estimates. 
Numerically, the gap can be used to control the rate of convergence of numerical 
computation methods such as discretization or finite element method by which one 
uses matrices to approximate a differential operator. The ability to solve for the 
first eigenvalue and eigenvector of these matrices is controlled by the size of the 
gap between the first eigenvalue and the rest of the spectrum. Understanding the 
behavior of the gap for collapsing convex polygonal domains is also relevant to 
computer graphics image rendering 
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In [25], using gradient estimates, the lower bound 7r^/4 was proven for the gap 
function. In [55], further generaHzations of the above result to Schrodinger opera- 
tors were obtained. More importantly, in that paper, the case when the potential 
function is not convex was discussed and the lower bound of the gap was obtained. 
After the Hessian of the log of the first eigenfunction was estimated, the result is 
the counterpart of the Li-Yau estimate [16^ on the first eigenvalue with Ricci cur- 
vature being bounded below. Note that studying the Schrodinger operators rather 
than the Laplacian is not only a generalization but is also necessary. For example, a 
simple proof of the log concavity of the first eigenfunction (A theorem of Brascamp- 
Lieb) was obtained in 25j by viewing the Laplacian as the limiting operator of a 
series of Schrodinger operators. The gap estimate was sharpened in [29.^ and [3T] : 
currently the best lower bound for the gap function on the Dirichlet Laplacian on 
a convex domain in M" is tt^ . A sharp upper bound for the gap function on convex 
domains in M" is given in Proposition 2 of [IT] . 

Our first result is a compactness theorem for the fundamental gap on triangular 
domains. Since ^ is invariant under scaling of the domain, we consider ^ on the 
moduli space of triangles which consists of all similarity classes of triangles. In terms 
of the smallest two angles, the moduli space of triangles may be represented by a 
triangle with one side removed in the a x (3 plane, where a point (a, /?) represents 
the similarity class of triangles with angles an < /Jtt < tt — an — /3tt; see Figure [T] 

Theorem 1. Let P he the set of all similarity classes of triangles. Then, ^ : P ^M. 



as defined in (0.1) is a proper map: 

Pc:^{TeP \ e(r) < C} 

is a compact set. 

This theorem shows that there exists a minimizer for the fundamental gap on 
triangular domains, and in particular, we expect the following. 



Conjecture 1. Let : P ^ R be defined by (0.1). Then, ^ has no extrema on 
the interior of the moduli space of triangles, and moreover, ^ is monotonic on the 
non- degenerate boundaries of the moduli space of triangles. 

The non-degenerate boundaries of the moduli space of triangles consist of isosce- 
les triangles, and the equilateral triangle lies at the non-degenerate vertex. Con- 
jecture 1 states that ^ has neither maxima nor minima for scalene non-degenerate 
triangles and moreover, for isosceles triangles with angles an = an < n — 2an or 
angles an < n/2 — an/2 = n/2 — an/2, ^{a) is strictly decreasing on a g (0, 1/3). 
The eigenvalues and eigenfunctions of the equilateral triangle T are explicitly com- 
putable [T3], and 

Consequently, Conjecture 1 would immediately imply the following conjecture, 
which also appeared in [2]. 

Conjecture 2. Let T C be a triangular domain. Then 
where equality holds iff T is equilateral. 



THE FUNDAMENTAL GAP CONJECTURE FOR POLYGONS 



3 



Our second theorem is a technical generahzation of the compactness theorem 
which provides conditions under which the fundamental gap remains bounded or 
becomes infinite as convex polygonal domains collapse to the unit interval. 

Theorem 2. Let {Qn}nen ^6 convex m-gons so that lim„_^oo Qn = [Oj !]■ Assume 
Qn C {{x,y) ^M? : x,y > 0}, and that the longest side Sn ofQn lies on the x-axis. 
Let the height hn of Qn be defined by 

(1.1) hn ■.= mf{b:Qn C [0,a] x [0,6]}. 

(1) // there exist rectangles i?„ D Q„ D r„ so that 

Area{Rn) - ^rea(r„) < 0(/ifJ, 

then ^(Qn) is bounded as n oo. 

(2) // there exists a convex inscribed polygon Un C Qn for which the following 
conditions are satisfied, then S,{Qn) —^oo as n oo. 

(a) The diameter of Un ^ as n —> oo. 

(b) One side S„ o/ [/„ is contained in Sn. 

(c) The height of C/„ = hn . 

(d) The height of Vn '■— Qn ~ Um satisfies h(Vn) < hn — 0{hn) for some 

In section 2, we study the behavior of the fundamental gap on triangular domains 
and prove the compactness theorem. In section 3, we provide examples of collapsing 
quadrilateral domains and prove our theorem for polygonal domains. Concluding 
remarks comprise section 4. 

2. The fundamental gap on triangular domains 

We consider here the Euclidean Laplacian on , which in rectangular and polar 
coordinates is respectively, 

A = ~{dl + dl), A = ~df. - r-^dr - r~^dl 

The Laplacian on a polygonal domain in M'^ with Dirichlet boundary condition has 
discrete spectrum 

< Ai < A2 < A3 < . . . . 

The vertex of a circular sector with the appropriate angle is a good local model 
for the corresponding vertex of a triangle [5j. This motivates one to approximate 
the eigenvalues of triangles using the variational principle and the eigenvalues for a 
circular sector, which can be explicitly computed using separation of variables for 
the Laplacian in polar coordinates. 

2.1. Eigenvalues for circular sectors. The eigenvalues of a circular sector of 
opening angle an and radius 1, with Dirichlet boundary condition are 

(2.1) {Afc,s}/c,seN = {{jk/a.s)^}k,seN, 

where ji,^s is the s*'' zero of the Bessel function of order v. Note that if the radius 
of the sector is r, the eigenvalues scale by r^^. 

We recall the variational or mini- max principle for the first eigenvalue. 



(2.2) Ai = inf 



\Vf\^dxdy 



/GDom(A):/7tO f^dxdy 
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where Dom(A) is the domain of A, so the infinium is taken over functions satis- 
fying the boundary condition. This is the infinium of the Rayleigh-Ritz quotient. 
Similarly, the second eigenvalue is 

(2.3) A2 = mf \ „, , — , 

/eDom(A):/#oj_L/i f^dxdy 

where fi is the eigenfunction for Ai, and orthogonality is with respect to These 
formulae are in [1| and A further property of the eigenvalues proven in |3] and 
[6] is domain monotonicity: if 17 C il' then for the Dirichlet eigenvalues 

Xkin) > Xk{n'). 



2.1.1. Asymptotic formulae for zeros of Bessel functions and eigenvalues of circular 
sectors. The following formulae are due to [14j, [18j, and [24J for Bessel functions of 
real order and date back to the work of [28] and [21] for Bessel functions of integer 
order. The first and second zeros of the Bessel function of order v are 

(2.4) ^•^_^^^__^^i/3 + o(^-i/3)^ , = 1^2, 

where ai is the i*"^ zero of the Airy function of the first kind so that 

(2.5) ai « -2.33811, and 02 w -4.08795. 

Consequently, we have the following estimates for the first two Dirichlet eigenvalues 
of the circular sector of opening angle air and radius one 

(2.6) \{a^)^\ + ^ + 0{a-^), z = l,2, 
where 

(2.7) ci = -ai22/3 w 3.71151827 and C2 -022^^^ w 6.48921613. 
We will also use the constaniQ 

(2.8) c'l = -a'i22/^ « 1.61722832. 



2.2. Gap behavior approaching the degenerate boundary of the moduli 
space of triangles. Since the gap function is invariant under scaling, we restrict 
to triangles with diameter one. Such a triangle has angles 

< aTT < /Jtt < Ti ~ ai: — /Jtt. 

The moduli space of all such triangles is itself a triangle in the a x f3 plane; see 
Figure [T] We are interested in the behavior of ^ on this moduli space, and in 
particular the behavior of ^ approaching the degenerate boundary of P which is 
the dashed vertical segment in Figure [l] 



^This constant arises from the asymptotic formula for the first zero of the derivative of the 
Bessel function which is related to the first Dirichlet eigenvalue of an obtuse isosceles triangle; see 

Eol. 
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(0, 1/2) 




(0, 0) 



Figure 1. Moduli space of triangles. 

2.2.1. Gap behavior away from the origin in the moduli space of triangles. Away 
from the origin in the moduli space of triangles, we are able to prove that ^ ~ a^^^^ 
as a ^ 0. Consider the triangle T with angles < air < (3tt < tt — an — Ptt, and 
assume for some fixed e > 0, /? > e. Let the side opposite an have length A, the 
side opposite /Jtt have length B, and the third side have length one. The law of 
sines states that 

sin(a7r) sin(/97r) sin(7r — an — Pn) 
A ^ B ^ 1 ■ 

Note that 

B= . f"(^") , ^1-Q(a)asa^0. 
sin(Q;7r + p7r) 

For a small, we may approximate the eigenvalues using two sectors, a larger sector 
of radius 1, and a smaller sector of radius 1 — 0{a), both with opening angle an; 
see Figure |2] By domain monotonicity, 

A2(T) > A2 (large sector) = X2{S), 

Ai(r) < Ai (small sector) = Ai(s). 

Then, 

m > X2{S) - Ai(.) ^ + - a^il-Oia)r " a^/^il-Oia^ + 

This shows that £,{T) behaves like a^'^^'^ as a 0. 

For some specific trajectories approaching the origin in the moduli space, we are 
also able to show that ^ ~ a~'^^^ as a ^ 0. If a, /3 ^ so that a = o(/3^), the 
following estimates demonstrate that ^ ^ a^"'/'^ as a — > 0. Consider the triangle 
with smallest angle an and vertices A,B,E. When x = \DE\ is very small we 
approximate f (T) using the larger triangle PBE, which is similar to the triangle 
ABD; see Figure [3j By domain monotonicity, 

X2iABE) > \2{PBE), and \i{ABE) < Xi{ABD), 
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Figure 2. Triangle with one collapsing angle and approximating 
circular sectors. 




Figure 3. Collapsing triangle and approximating triangle. 

so that 

X2iABE) - Xi{ABE) > X2{PBE) - Xi{ABD) = X2iABD){l + x)"^ - XiiABD). 
Since 

X^iABD) =\ + ^ + O(a-i), and XMBD) = \ + ^ + 0{a-'), 
for X < 1 we expand 

For 

v2/3 



C2Q! 

X < — — 



we then have the estimate 



X2{ABE) - Xi{ABE) > 



,4/3' 

for some ^ > 0. 

2.2.2. Gap behavior approaching the origin along a non- degenerate boundary in the 
moduli space. If we approach the origin in the moduli space along trajectories near 
the non-degenerate boundary which meets the origin, these triangles are "almost" 
isosceles, and we can show that ^ ~ a~'^^^ as a — > 0. First, we consider obtuse 
isosceles triangles with diameter one. By [10 Theorem 5.3, the first eigenvalue is 

The second eigenvalue of T was also computed in [10] , 

A2(T)=:i + ^ + 0(a-^/3). 



THE FUNDAMENTAL GAP CONJECTURE FOR POLYGONS 



7 



Then, 

which shows that S,{T) ~ a"^/'^ as a — + 0. Note that [10 computed further terms 
in the asymptotic expansion of ^ in this case. 

When two angles are collapsing at approximately the same rate so that the 
triangle is "almost" isosceles, we estimate using the variational principle. Let the 
two collapsing angles be air and (3Tr. Assume that there exists a constant c € (0, 1), 
to be specified at the end of this argument, so that 

Q. 

a < P, c < — < 1 as a, /3 — > 0. 

Let the third angle of the triangle be 7 with opposite (longest) side length 1, and 
let A and B be the sides opposite angles air and /Stt, respectively. In Figure |4] 
this is triangle PRT. Consider the isosceles triangle with angles /?7r, /Stt, tt — 2/37r, 
contained in the original triangle; in Figure |4] this is triangle QRT. The sides of 
this triangle are A,2Acos{j3Ti) . Let Ai(/37r) be the first Dirichlet eigenvalue for 
this triangle and let Ai be the first eigenvalue of the original triangle. By domain 
monotonicity, 

^ ^'^^^^ - (AcosilnWP^ + (Acos(/3t))^/^V3 + ^^^-'^'^ 

By domain monotonicity, we approximate A2 from below by the second eigenvalue 
of the isosceles triangle with angles a7r, aTr, tt — 2a7r and sides B, B,2B cos{aTr), 
which contains the original triangle and is triangle PST in Figure [4j The second 
eigenvalue for this isosceles triangle as computed by [10] is 

A2(a7r) = — / „ + 7^ . ^\x2 4/. + 0(a"^/^). 

(i3cos(a7r))^a^ (B cos(a7r))"'a^/'^ 

Therefore, 



(Bcos(a7r))2a2 (B cos (aTr)) 2^4/3 (^cos(/37r))2/32 (Acos(/37r))2/34/3 

By the law of sines, 

Asin(/37r) = Bsin(a7r) ^ A^^ Ba + 0{(3^), 

when a and (3 are small. Moreover, cos(a7r) = 1 + 0{a^), and cos(/37r) = 1 + 0{(3^) 
as a, (3 ^ 0. This shows that we need to approximate 

Cl Cl 

^2^4/3 A2/?4/3' 

By the hypothesis, 

c/3 < a < /3, 

so that when a and /3 are small, the law of sines gives 

cB <A<B. 

Then, 
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T 




Figure 4. Triangle with two angles collapsing and approximating 
isosceles triangles. 




B p., P1 D P2 P2' E C 



Figure 5. Arbitrary collapsing triangle. 

So, we require 

Since c < ^ < 1 , this will be satisfied for 

c > i-^j « 0.660100001. 

This shows that for collapsing triangles with c < ^, ^(T) is bounded below by a 

constant multiple of a"^^"^. 

Due to these calculations and the work of we expect that ^ is a polyhomo- 
geneous function on the moduli space of triangles blown up at the origin. Such a 
regularity result for the gap function is an interesting open problem for polygonal 
domains, and in particular, it would show that for any family of collapsing triangles 
with smallest angle a, ^ oT^I"^ as the triangles collapse. 

2.3. Proof of the compactness theorem. Assume the smallest angle of the 
triangle AB = aTr; see Figure [s] We wish to estimate \2(^ABE) — Xi{ABE) from 
below. Assume \DE\ < \BD\ = 1, and fix 

2 

< e < -. 
9 

Let Pi be between B and D so that \PiD\ = a'^, and P[ be between B and Pi so 
that |P{D| = 2a^ If \DE\ > 2a\ let P2 be between D and E so that \P2D\ = a\ 
and P2 be between P2 and E so that \P2D\ = 2a^ If \DE\ < 2a% we do not define 
or use the points P2,P2- If \DE\ > 2a^ let U be the trapezoid AQ1P1P2Q2 and 
similarly let [/' = AQ[P{P^Q'^. Let V = ABE ~ U and let V = ABE - U'; see 
Figure |t] If IL'^I < 2a% we let U = AQiP^E, U' = AQ[P{E, V = ABE - U 
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and V' = ABE — U' . In the estimates to follow, we show that we may estimate 
X2iABE) - Xi{ABE) using A2(?7') - Ai(C/')- 

Let fi be the eigenfunction for Xi{ABE), i ~ 1,2. The height of V is at 

most 

(1 — a'^) tan(Q;7r) sa (1 — a'^)7ra. 
By the one dimensional Poincarc Inequality 

Ivl^M'^ lul^M' ^ andi^^^^>^. 



Assume 

/ /2 = 1, and / = (3. 
Jabe Jv 
By the variational principle. 

Therefore, 

For simplicity in the arguments to follow, we replace all constant factors multiplying 
positive powers of a by a constant factor of 1, since no generality is lost as a — *■ 0. 
Let p he a, smooth compactly supported function so that 

(2.9) P\u = h P\v'=0- 
We may choose p so that 

(2.10) |Vp|<^, and |Ap|, |A(p2)| < 

For the arguments to follow, we use the sign convention for the Euclidean Laplacian 
so that —A has positive spectrum. Note that 

(2.11) - (p/»)A(p/,) - Kp^ff - ffpAp - 2/,p(Vp)(V/0. 

2.3.1. Estimate for Xi(U'). Since p vanishes on the boundary of [/', pfi is an ad- 
missible test function for the Rayleigh quotient on U' (2.2) which we may use to 
estimate Xi{U') from above. By ( |2.11| , 



Ai(f/')< Ai + 



{-pApff - 2pVp/iVA) 



Since 

and Vp, Ap = on U, we have 

*-2 < JL < Q,2/3-3e 



/ {-pApf!~2pWphWh)<^ f fl<, <,a- 
Ju' ^ JU'-U 

Noting that 

Ju' Ju 

we then have 

(2.12) Ai([/0 < Ai+ " ,,3 ^ < Al + a^/3-^^ 

1 — a'^i'^^'^ 
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This gives the required estimate for the first eigenvalue. 



2.3.2. Estimate for A2 (?/'). Since p/2 is not a priori orthogonal to the first eigen- 
function for U' , we must modify it to use the Rayleigh quotient (2.3 1 to estimate 
X2{U'). Since pfi is not orthogonal to the first eigenfunction for U' because both 
are positive, there is some a S M such that p/2 + apfi is orthogonal to the first 
eigenfunction for U' . We may then use p/2 + ap/i as a test function for the Rayleigh 
quotient on U' . Integrating by parts, 

(2.13) / |Vp/,p = -/ pfAipf.) = Kf p'f!~l [ ^p'^f--[ P^Pfl 

JU' JU' JU' ^ JU' JU' 

We estimate 



iv(p/.: 



p^ft 



< 



\p^p\fi 



(2.14) 



< a" 



f2 ^ -2c 



U'-U 



2/3-3e 



since Ap and Vp vanish identically on U. We compute that J^^, V(/o/i)V(p/2) = 

- / pfiAiph) = \2 I pVi/2 - 2 / pVph'^f2 - f p/lAp/2 
JU' JU' JU' JU' 



(2.15) 

and/^,V(p/i)V(p/2 
(2.16) 

This gives the inequality 
2 



p/2A(p/i) = Ai /^^ p'hh - 2 / PVP/2V/1 



U' 



P/2AP/1 



V(p/i)V(p/2) - (Ai + A2) / pVi/2 

U' 



< 



U'-U 



1/1/2 

v2e 



pVpV(/i/2; 



U'-U 



(2.17) < (^1^ I/1/2I + 2 1^ I/1/2I) < 

which follows from integration by parts and the Schwarz inequality. Expanding, 
/ |V(p/2 + ap/i)|2-A2 / (p/2 + ap/i)' -/ + // + ///, 

Ju' JU' 



where 



and 



|Vp/2 



|Vp/i 



p'f? 



111= f 2a(Vp/i)(Vp/2)-A2 / 2apVi/2. 

JU' JU' 

By ( [2l4t , 

(2.18) / < a2/3-3.^ 
To estimate II, we note that Ai < A2 so that 

(2.19) II<a'(^j^ \Vph\^-X,j^ p'fl^ <aV/3-3% 
which follows from (2.141. 
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Note that by the orthogonahty of /i and /2 and the Schwarz inequahty, 



(2.20) 



p'/i/2 



ABE 



(1 - p')/l/2 



< 



/1/2 



< a 



2/3-e 



By (2.171, (2.201, and adding and subtracting Ai Jjj, ap^/1/2, we estimate III, 



/// < \a\ 
(2.21) 



(Vp/i)(V/7/2)-(Ai + A2) / p^hf: 



+ |a|(A2-Ai) 



P''/l/2 



< |a|(a2/3-3^ + (A2-Ai)a2/3-). 



By (2.201, and since J^, /f > 1 - a 



2/3- 



(2.22) 



(p/2 + aphf > (1 + a'){l - a^/^-^) - 2\a\a'/^-'^. 



We now estimate the Rayleigh quotient for p/2 + apfi using (2.181, (2.191, (2.21 1, 
and (|2.22|) 



A2(f/') < A2 + 
which shows that 



(1 + a2)(l - a2/3-e) _ 2|a|a2/3-3e 



A2(C/') < A2 + + (A2 - Al)a2/3-^ 

2.3.3. Gap estimate. Using our estimates for Ai([/'), 

A2 - Ai > A2(C/') - Ai(C/') - a2/3-3e _^ _ Ai)a2/3-e _ ^2/3-3e 

We then have 

(A2 - Ai) > ^3^57337 (^2(f/') - Ai(C/')) - a^/^-^S 
which shows that 

A2 - Ai > (A2(C/') - Ai(t/')) - 0(a2/3-3.)_ 
By the main theorem of ^25j, since the diameter of U' is at most 4a', 



A2(C/')-Ai(C/')> 



which shows that 

A2 - Ai > Ca-^' 

and is therefore unbounded as a ^ 0. We have shown that for any triangle with 
one or two small angles, ^(T) becomes unbounded as the small angles collapse. 
Therefore, any sequence of triangles for which f (T) is bounded above by a constant 
C cannot contain a collapsing subsequence. Since it must contain a subsequence 
of triangles T„ for which ^(T„) — > C, it must therefore contain a subsequence 
converging to a triangle T for which ^(T) = C. This completes the proof of the 
compactness theorem. □ 
Remark: The numerical methods of [20] use linear combinations of Bessel 
functions to accurately approximate the Dirichlet eigenvalues of polygonal domains. 
One may use such methods together with our compactness theorem to numerically 
prove a lower bound for the gap function on all triangular domains. 
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3. Fundamental gap on polygonal domains 

Some readers may find our compactness result for triangles counterintuitive be- 
cause of the familiar example of collapsing rectangles. 

3.1. The fundamental gap on quadrilateral domains. For a rectangle, R, 
one may compute the Dirichlet or Neumann eigenvalues and eigenfunctions using 
separation of variables. For a rectangle with side lengths a, b where we assume 
h < a, the first two Dirichlet eigenvlaues are 

Al = TT^ 

The fundamental gap is then 















^2 = tM^- 





The gap conjecture becomes sharp when the rectangle collapses as 6 ^ 0, and we 
note that the fundamental gap for the interval [0, a] is 

C([0,a])=37r^ 

This is intuitive since rectangles collapse uniformly to the segment; consider the 
domain the fixed square [0,1] x [0,1] with coordinates {x,t) where t = by. The 
Laplacian in coordinates {x,t) is related to the Laplacian on the rectangle [0, 1] x 
[0,6] with coordinates {x,y) by 

^dl-dl = -dl-b'dl 

As 6 ^ 0, the operator converges in some sense to the operator —9^ on the interval, 
and the gap converges to the gap on the interval [0, 1]. 




Figure 6. Collapsing quadrilateral with bounded gap. 



For triangles, the collapse is not uniform; points opposite the longest side col- 
lapse at diS^erent rates. Points near the vertex opposite the longest side collapse 
"more slowly" in some sense than points near the smallest angle. For general polyg- 
onal domains, there is a subtle relationship between uniformity of collapse and the 
behavior of the gap: the examples in Figure [6] and Figure [7] are both approximately 
rectangular yet in the first example the gap remains bounded as the quadrilateral 
collapses while in the second example the gap becomes unbounded. This follows 
from our second theorem. 
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Figure 7. Collapsing quadrilateral with unbounded gap. 

3.2. Proof of Theorem 2: bounded gap. The first case follows almost imme- 
diately from domain monotonicity estimates. For ease of notation, we drop the 
subscript n. The rectangle R = [0,A] x [0,B]. Since Q is collapsing to [0,1], we 
must have B k h and A ^ \. The rectangle r = [0,a] x [0,6]. By the variational 
principle, 

e(Q)< A2(r)-Ai(i?) = 7r2 ^ ^ ^ 



Since the area of ~ /i, by the hypothesis the area of r « ft, + 0{h^). Since a < A, 
and b < B, it follows that a — > 1. We calculate 

The last two terms are bounded since a,^ — > 1. Since 

B^-b^ hO{h^) 
{bBy W 
is bounded as /i — > 0, we see that ^ is bounded as n — > oo. 




Figure 8. Schematic diagram of collapsing polygon with un- 
bounded gap. 

3.3. Proof of Theorem 2: unbounded gap. We generalize our estimates for 

arbitrary collapsing triangles. Let ?7' D f/ be a convex inscribed polygon so that one 
side r of U' satisfies 5 D T D S and |r - S| = ft'', where 5 G (0, 1) will be specified 
later. We can define such an inscribed polygon since the diameter of ?7 — > 0, and the 
length of the longest side of Q ^ 1 . Note that with these hypotheses the diameter 
of i7' ^ as ft ^ 0. Let /j be the eigenfunction for Q with eigenvalue Aj, for 
i = 1,2. By convexity, Q contains an inscribed right triangle T of height ft = ht{U) 
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and base at least |S'|/2, where \S\ is the length of the longest side of Q. Scale Q 
so that the base of T is one; no generality is lost in showing ^ is unbounded. For 
illustration, in Figure [s] Q = AEFJL, T ^ JCE, h = \ JC\, and U = JKBDH. 
Assume /j are normalized so that 



= and let / = /3. 



By the one dimensional Poincarc inequality, 



> 



-{ht{v)Y^ 

and 

since U and U' have the same height. For w 0, the measure of the smallest angle 
of T is approximately h = ht{U). By domain monotonicity, 

Estimating (i as we did for triangles, 
/37r2 7r2(l-/3) 



< 



/jV/.P+/|V/.P = A. = A.(T)<^ 



{ht{v)Y' {ht{u)Y -Jy' 'Ju' ' - {ht{uw ' {ht{u)Y/^- 

This gives the following estimate for /3, 

ht{Uf/^ht{VY 



(3.1) /3< 



We have dropped the constant factor. By considering the leading asymptotic be- 
havior as h ^ Q, we have the estimate 

(3.2) (3<h^l^-''. 
Let 

(3.3) e=\-x>Q. 

3.3.1. Estimates for Ai. We prove an estimate of the form 

Ai(C/') < Al(Q) + /^^ 



for some fixed y > 0. Define the cut-off function p as in (2.9 1 so that, 

|Vp| < h-\ |Ap| < h'"^^. 



Using the test function pfi in the Rayleigh quotient for U' , by (2.131, 

(3.4) Ai(C/') < Ai(Q) + < Ai(Q) + h^-^' . 

We again absorb all constant factors multiplying positive powers of h into a constant 
factor of 1, since this does not change our limiting estimates for h ^ Q. 
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3.3.2. Estimates for A2. Wc now require an estimate of the form 

X^iU') < A2(Q) + hy + (A2(Q) - XiiQ))hy' 

for some y,y' > 0. Since the function pf2 is not a priori admissible as a Rayleigh 
quotient test function for A2(C/') we again modify it to make it orthogonal to the 
first eigenfunction on U' . So, we consider the test function 

/0/2 + apfi. 



Using (2.131, we estimate 



|V(p/,)p-A, / 



/ pVf < I I Ap'f^ + f \pAp\ff 
Ju' ^ Jv Jv 



(3.5) 



< h 



-25 



-2S 



U'-U 



-25 



V 



since Ap and Vp vanish identically on U. By (|2.15|) and (|2.16|, 
2/ V(p/i)V(p/2) - (Ai + A2) / 



U' 



U' 



< 



1/1/2 



U'-U 



pVpV(/i/2) 



U'-U 



(3.6) 



< h- 



1^1/1/21 + 2 /^l/l/2|) <h^-'', 

which follows from integration by parts and the Schwarz inequality. Expanding, 
/ |V(p/2 + ap/i)|2-A2 / (ph + aph)^ = 1 + 11 + 111, 

Ju' JU' 



where 



and 



|Vp/2 



|Vp/j 



111= I 2a(Vp/i)(Vp/2)-A2 / 2apVi/2. 

JU' JU' 



By our calculations for triangles, our estimate for /3, and the estimates (3.5) and 
(Pl, 



(3.7) 
(3.8) 
(3.9) 

Moreover, 
(3.10) 



// < a^h"-^^, and 
/// < \aW-^' + (A2(Q) - X,{QW-^'). 

(p/2 + aphf > (1 + a^){l - h') ~ 2\a\h 



■-2S 



Using these estimates, we estimate the Rayleigh quotient for p/2 + apfi, 

(3.11) A2(C/)<A2(Q) + (i + „2)(i_/,.)_2|a|/,e-25 ■ 

By considering the behavior as h ^ 0, 

X2{U') < A2(Q) + h'-^^ + h'-^\X2iQ) - Ai(Q)). 
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3.3.3. Gap estimate. Using our estimates for Xi(U'), 

HQ) - Ai(Q) > A2(C/') - Ai([/') - h'-^' + (A2(Q) - Xi{Q))h'-^', 
which shows that 

(3.12) HQ) - HQ) > ''^"".'ll'lP - h^-''- 

Since we are free to choose any S > 0, we may choose 5 e (0, e/2). The diameter of 



[/' is therefore vanishing as 0, so (3.121 and the main theorem of imply 



HQ) ~ Ai((3) — > cx), as /i ^ 0. 

□ 

Remark: This technical theorem shows that the gap function is sensitive to 



the rate at which boundary points converge to the longest side. If the height ( 1.1 ) 
of the polygon is /i, and if all boundary points are collapsing at the same rate with 
an error controlled by 0(h?)^ then the gap remains bounded. However, if the rate of 
collapse varies by 0{h^) for some a; < 5/3, then the gap becomes unbounded. For 
5/3 < a; < 3, our methods do not give information on the gap behavior, however for 
a generic polygonal domain, we expect that the rate of collapse varies by at least 
0{h^/^^'). One may conclude from Theorem 2 that generic polygonal domains 
which collapse to the interval have unbounded gap. This is consistent with [3] and 
^12, : in particular, see the remarks after Corollary 2.1 of 



4. Concluding remarks 

We have analyzed the behavior of the gap on the degenerate boundary of the 
moduli space of triangles and generalized our methods to polygonal domains. Sup- 
ported by numerical estimates, we conjecture that the gap minimizer for triangular 
domains is the equilateral triangle. Further analysis of the behavior of ^ on the 
interior of the moduli space and its non-degenerate boundaries is a reasonable ap- 
proach to prove the gap conjecture for triangular domains; this is an open problem. 
For convex polygonal domains of more than three sides, our examples indicate that 
generically one expects the gap to become unbounded in the moduli space of convex 
m-gons approaching the degenerate boundaries. Some reasonable open questions 
are the following. Is it possible to determine necessary and sufficient conditions 
under which the gap remains bounded as convex polygonal domains collapse to the 
interval? Among non-collapsing domains, what is the gap minimizer, or does it 
exist? On the moduli space of convex m-gons, may the gap function have interior 
extrema, and what are its regularity properties? We hope to inspire readers to 
attack some of the many open spectral problems, and that our work is a useful 
contribution to understanding the fundamental gap. 



Acknowledgements 

The first author is partially supported by the NSF career award DMS 0347033. 
The second author would like to thank Vincent Grandjean, Rafe Mazzeo, Richard 
Melrose, and Michael Taylor for helpful discussions and encouragement, and the 
Mathematical Sciences Research Institute for a stimulating and beautiful environ- 
ment in which a portion of this work was completed. 



THE FUNDAMENTAL GAP CONJECTURE FOR POLYGONS 



17 



References 

[I] M. Abramowitz and I. Stegun, Handbook of mathematical functions with formulas, graphs, 
and mathematical tables, National Bureau of Standards Applied Mathematics, Series 55, for 
sale by the Superintendent of Documents, U.S. Government Printing OfRce, Washington, D.C., 
1964. 

[2] P. Antuncs and P. Frcitas, A numerical study of the spectral gap, J. Phys. A: Math. Theor. 

41, 2008, 055201, 19 pages. 
[3] D. Borisov and P. Freitas, Singular asymptotic expansions for Dirichlet eigenvalues and eigen- 

functions on thin planar domains, to appear in Ann. Inst. H. Poincare Anal. Non Lineaire, 

arxiv math.SP 0712.3245vl. 
[4] I. Chavel, Eigenvalues in Riemannian Geometry, Academic Press, 1984. 

[5] J. Cheeger, Spectral Geometry of Singular Riemannian Spaces, J. Diff. Geo., 18, (1984) 575- 
657. 

[6] R. Courant and D. Hilbert, Methods of Mathematical Physics, Volume 1, Interscience Pub- 
lishers, 1937. 

[7] S. Finch, Airy Function Zeroes, supplementary material for Mathematical Constants, Cam- 
bridge University Press, 2003. 

[8] S. Finch, Bessel Function Zeroes, supplementary material for Mathematical Constants, Cam- 
bridge University Press, 2003. 

[9] J. D. Foley, A. van Dam, S. K. Feiner, J. F. Hughes, Computer Graphics: Principles and 
Practice in C, Addison- Wesley Professional, 1995. 

[10] P. Freitas, Precise bounds and asymptotics for the first Dirichlet eigenvalue of triangles and 
rhombi, J. F. A. 251, (2007), 376-398. 

[II] P. Freitas and D. Krejcirik, A sharp upper bound for the first Dirichlet eigenvalue and the 
growth of the isoperimetric constant of convex domains, Proc. Amer. Math. Soc. 136, 2008, 
2997-3006. 

[12] L. Friedlander and M. Solomyak, On the Spectrum of the Dirichlet Laplacian in a Narrow 

Strip, preprint arxiv math.SP 0705.4058vl. 
[13] J. R. Kuttler and V. G. Sigillito, Eigenvalues of the Laplacian in two Dimensions, SIAM 

Review, 26, no. 2, (1984) 16.3-193. 
[14] T. Lang and R,. Wong, "Best possible" upper bounds for the first two positive zeros of the 

Bessel function .J,^(x): the infinite case, J. Comp. Appl. Math., 71, 1996, 311-329. 
[15] N. Lebedcv, Special functions and their applications, Prentice-Hall, Inc., 1965. 
[16] P. Li and S. T. Yau, Estimates of eigenvalues of a compact Riemannian manifold. Geometry 

of the Laplace operator (Proc. Sympos. Pure Math., Univ. Hawaii, Honolulu, Hawaii, 1979), 

pp. 205-239, Proc. Sympos. Pure Math., XXXVI, Amer. Math. Soc, Providence, R.I., 1980. 
[17] R. Liboff, Nodal and other properties of the second eigenfunction of the Laplacian in the 

plane. Quarterly of Applied Mathematics, Vol. LXIII, No. 4, 2005, 673-679. 
[18] L. Lorch and R. Uberti, "Best possible" upper bounds for the first two positive zeros of the 

Bessel functions - the finite part, J. Comp. Appl. Math., 75, 1996, 249-258. 
[19] A. Mclas, On the nodal line of the second eigenfunction of the Laplacian in R^, J. Diff. 

Geometry, 35, 1992, 255-263. 
[20] C. Molcr, Accurate bounds for the eigenvalues of the Laplacian and applications to rhombical 

domains, Technical Report No. CS 121, Computer Science Department, Stanford University, 

1969. 

[21] F. W. J. Giver, Bessel Functions Part II, Zeros and Associated Values, Royal Society Math 

Tables Vol. 7, Cambridge University Press, 1960. 
[22] G. Polya, On the characteristic frequencies of a symmetric membrane. Math. Z. 63, 1955, 

331-337. 

[23] G. Polya and G. Szego, Isoperimetric inequalities in mathematical physics, Ann. Math. Stud- 
ies, 27, Princeton Univ. Press, 1951. 

[24] C. K. Qu and R. Wbng, "Best Possible" Upper and Lower Bounds for the Zeros of the Bessel 
Function Jv{x), Transactions of the American Mathematical Society, Vol. 351, No. 7, 1999, 
2833-2859 

[25] I. M. Singer, B. Wong, S.T. Yau, and S. S. T. Yau, An estimate of the gap of the first two 
eigenvalues in the Schrodinger operator, Ann. Scuola Norm. Sup. Pisa Gl. Sci. 4, Vol. 12, No. 
2, 1985, 319-333. 



18 



ZHIQIN LU AND JULIE ROWLETT 



[26] B. Siudeja, Sharp Bounds for Eigenvalues of Triangles, Michigan Matli J., Vol, 55, Issue 2, 
2007, 243-254. 

[27] M. van dc Berg, On the condensation of the free Boson gas and the spectrum of the Laplacian, 

J. Stat. Phys., 31, 1983, 623-637. 
[28] G. N. Watson, A Treatise on the Theory of Bessel Functions, Cambridge University Press, 

1966. 

[29] S. T. Yau, An estimate of the gap of the first two eigenvalues in the Schrodinger operator. 
Lectures on Partial Differential Equations, 223-235, New Stud. Adv. Math., 2, Int. Press, 
Somcrville, MA, 2003. 

[30] S. T. Yau, Nonlinear Analysis in Geometry, Enseignement Math., 33, 1987, 109—158. 
[31] Q. H. Yu and J. Q. Zhong, Lower bounds of the gap between the first and second eigenvalues 
of the Schrdinger operator, Trans. Amcr. Math. Soc. 294, no.l, 1986, 341-349. 

Department of Mathematics, 103 MSTB - University of California, Irvine, CA 92697- 
3875. 

E-mail add/i'ess: zlu@math.uci.edu 

Department of Mathematics, South Hall 6607, University of California, Santa Bar- 
bara, CA 93106. 

E-mail address: rowlett9math.ucsb.edu 



